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Abstract 

We estimate the relative probabilities of ^'(2S) to J/^f production at BNL-RHIC 
and T(nS) production at the LHC and Fermilab in p-p collisions, using our recent 
theory of mixed heavy quark hybrids, in which the ^'(2S) and T(3S) mesons have 
approximately equal normal qq and hybrid qqg components. 

PACS Indices:12.38.Aw,13.60.Le,14.40.Lb,14.40Nd 

1 Introduction 

There has been a great deal of interest in the production and polarization of heavy quark 
states in proton-proton collisions. This was motivated in part by the J/ty,^' production 
anomalyfl], in which the charmonium production rate was larger than predicted for J/^f, 
and much larger for than theoretical predictions in proton-proton (p-p) collisions; and 
the disagreement between experimental measurements of polarization of J/^f,^>' [2] states 
produced at high energy and theoretical predictions [3]. In addition to being an important 
study of QCD, it also could provide the basis for testing the production of Quark Gluon 
Plasma (QGP) in relativistic heavy ion collisions (RHIC). We use the notation E=i/i. 

At the proton-proton (p-p) energies of the Fermilab, BNL-RHIC, or the Large Hadron 
Collider the color octet model[3, 4, 5] dominates the color singlet model, as was shown in 
studies of J/^> production at E=200 GeV at BNL[6, 7]. In this color octet model heavy quark 
state production cross sectons require partonic distributions and nonperturbative matrix 
elements. In our present work we use this framework, but some of our nonperturbative 
matrix elements are taken from a recent publication [8] in which it was shown that some of 
the charmonium and upsilon states are mixed meson and hybrid meson states, which followed 
the work in which it was shown that these states are not pure hybrids [9]. This is particularly 
important for the application of the octet model, since the heavy quark QQ in a hybrid are in 
a color octet representation. The mixed hybrid model of Ref[8] can also explain the ty'(2S) 
production anomaly, with a much larger production of ty'(2S) in high energy collisions than 
standard model predictions[10], the famous p — n puzzle, and the anomolous production 
of sigmas in the decay of T(3S) to T(15')[ll]. In section 2 we discuss mixed heavy quark 
hybrids and the solution to these puzzles. 
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In the present work on heavy quark state production in p-p collisions we estimate the 
relative production of W(2S) to J/* at BNL-RHIC (200 GeV)[6, 7], and the relative pro- 
duction of (T(2S) + T(3S)), T(1S) states at the Large Hadron Collider (2.76 TeV)[12], and 
(T(nS) states at Fermilab (38.8 Gev)[13, 14]. 

We closely follow the formulation of Nayak and Smith[15] in which they calculated J/\I/ 
and ty'(2S) production, produced by unpolarized or polarized p-p collisions at yfs = 200 GeV 
for helicity A = and A = 1. We generalize this by considering the production of a heavy 
quark state $, where $ is Charmonium J/\E',\I/ / or Bottomonium T(nS). It is important 
to note that in Ref[8] the ty'(2S) and T(35) were found to be approximately 50% a heavy 
quark meson and 50% a hybrid. 

2 Review of mixed hybrid heavy quark mesons 

The Charmonium and Upsilon (nS) states which we are studying are shown in Figure 1 



" C 3 S > i < >, :^ fs fs 



1 0.0313 



J-' > 4'<:is> i « >< >~7 y<is> v>- >< > 



Figure 1: Lowest energy Charmonium and Upsilon states 



2.1 Heavy quark meson decay puzzles 

Note that the standard model of the ip'(2S) and T(3<S) as cc and bb mesons is not consistent 
with the following puzzles: 

1) The ratio of branching rarios for cc decays into hadrons (h) given by the ratios (the 
wave functions at the origin canceling) 

B(g'(cc) -» h) B(V(cc) -» e+e-) ^ 

5(J/>(cc) ->• /i) £(J/*(cc) ->• e+c) ' ' 

the famous 12% RULE. 

The p — 7r puzzle: The ^ / (2S f ) to J/\l/ ratios for p — n and other h decays are more than 
an order of magnitude too small. Many theorists have tried and failed to explain this puzzle. 
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2) The Sigma Decays of Upsilon States puzzle: The a is a broad 600 MeV ir—ir resonance. 

T(2s) — > T(1S) + 2n large branching ratio. No a 
T(3s) — > T(1S) + 2n large branching ratio to a 

We call this the Vogel An = 2 Rule [11]. Neither of these puzzles can be solved using 
standard QCD models. We have solved them using the mixed heavy hybrid 

2.2 Hybrid, mixed heavy quark hybrid mesons, and the puzzles 

First, we used the method of QCD Sum Rules, starting with the vector hybrid current 
(J PC = 1 ) for heavy quark hybrid mesons 

A™ 

T ^ a a b /-in b 

JhHv = q A l IS^^uvlB > 

where A,B are flavor indices, a,b are color indices, and G™^ is the gluon field operator. The 
QCD Sum Rule method starts with the two-point correlator for a heavy hybrid 

n^(g 2 ) = i J d'xe^ < 0\T[r HH (x)r HH (0)]\0 > . (1) 

Using the Jhh^ hybrid current and the QCD Sum Rule method it was shown[9] that no 
low-lying charmonium state is a hybrid. 

Then the mixed vector (J PC = 1 ) charmonium, hybrid charmonium current 

J" = bJ^ + Vl- b 2 J^^with 

J H = <? c 7 75? c , 

where is the current for a 1 charmonium state, and J^ HH mu is the heavy 
charmonium hybrid current given above, was tried [8]. A solution satisfying all 
QCD Sum Rule conditions was found for the ty'(2S) to be a mixed heavy hybrid, 
with b ~ —.7 

|tf'(2s) > = -0.7\cc(2S) > +Vl-O.5|cc0(2S') > . (2) 

Using this method for bottom quark mesons, it was found that the T(3S) 
state was a mixed hybrid: 

\T(3S) > = -0.7166(35) > +VT^0l\bbg{3S) > . (3) 

Therefore for both the W(2S) and T(3S) b=-.7, and both have a 50% proba- 
bility of being a hybrid and 50% probability of being a standard qq meson. In 
Ref[8] it was shown that this solves the ty'(2S) to J/ty the p — n puzzle, and the 
Vogel An = 2 Rule. 

We now investigate how this mixed heavy hybrid meson theory affects the 
production of charmonium and upsilon states via p-p collisions. 
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3 Unpolarized p-p collisions 

The cross sections for J/ty and ^>'(2S) production in the color octet model 
are based on the cross sections obtained from the matrix elements for quark- 
antiquark and gluon-gluon octet fusion to a hadron H, cr qq ^.H(\) and a gg ^H(\), with 
A the helicity, as illustrated in Figure 2. 




Figure 2: Gluon and quark- ant iquark color octet fusion producing hadron H 



In Ref[15] the matrix elements used were derived by Braaten and Chen[4]. 
The three octet matrix elements needed are < Of( 1 S , ) >, < Of( 3 Si) >, and 
< Of ( 3 P ) >, with $ either J/tf, V'(2S), or T(nS). Since these matrix elements 
are not well known, Nyyak and Smith[15] use three scenerios: 1)< Of ( 1 S , ) >=< 
Of ( 3 P ) > /m 2 =.0087, 2)< Of (%) >=.039 and < Of ( 3 P ) >=0, 3)< Of (%) >=0 
and < Of ( 3 P ) > /m 2 =. 01125, with < Of ( 3 Si) >=.0112 in all scenerios and all 
having units GeV 3 . Note that these matrix elements are not used to obtain the 
wave functions of the heavy quark meson states. 

To obtain the production cross sections one needs to multiply by the quark- 
antiquark or gluon parton distribution functions, giving 

o"pp^*(a) = / — f q (x,2m)f q (a/x,2m)a qq ^ H (\) 

J a X 

+f g (x, 2m)f g (a/x, 2m)a gg ^ HW , (4) 

where a = Am 2 / 's, with m — 1.5 GeV for charmonium, and 5 GeV for bottomonium. 
f g (x,2m), f q (x,2m) are the gluonic and quark distribution functions evaluated at 
Q = 2m. First we consider unpolarized p-p collisions, using[15] scenerio 2. With 
scenerio 2 matrix elements the production cross sections[4, 15] for $ for helicity 
A = and 1 are 



o- pp ^$ (A=0 ) = A® / — f g (x,2m)fg(a/x,2m) 

J a X 

o"pp^$(a=i) = A$ / — [f g (x,2m)fg(a/x,2m) + 0.613((/ d (x, 2m)/j(a/x, 2m) 

J a X 
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with A® = 288m3 s < Of^So) >, a = Am 2 /s; where m = 1.5 GeV for charmonium, 
and 5 GeV for bottomonium. f g (x,2m), f q (x,2m) are the gluonic and quark 
distribution functions evaluated at Q = 2m. 

The main purpose of the present work is to explore the effects of matrix 
elements for ty'(2S) and T(nS) with n=2,3. We compare our results with the 
hybrid model to the standard model. In the standard model the states are (nS) 
quark-antiquark states, and the ratios of the matrix elements for n greater than 
1 is given by the squares of the wave functions. Note that the basis for the 
octet model being used is the nonrelatavistic QCD model[3, 4, 5], with a model 
potential for the quark anti-quark interaction giving bound states. A harmonic 
oscillator potential can be used to approximately give the energies of the first 
few states, which is what is needed in the present work. For the octet matrix 
elements illustrated in Figure 1, however, one must use QCD directly, and we 
use the results of Refs.[3, 4, 5, 15] for these matrix elements. 

To approximate the ratios of matrix elements in a nonrelativistic quark model 
for these heavy quark meson states we use harmonic oscillator wave functions [16], 
with $(15) = 2Exp[-r/a ]/a 3 J 2 , $(25) = $(15)(1 - r/a )/2 3 / 2 , and $(35) = $(15)(1 - 
2r/3a + 2r 2 /27a 2 )/3 3 / 2 . Defining Nl= / |$(25)| 2 divided by / |$(15)| 2 for the 2S to 
IS probability, and simillarly N2 for the 3S to IS probability, we find Nl=0.039, 
N2=0.0064, N3=N2/N1=.16. This is a very rough estimate. The cross sections 
for the mixed hybrid states are enhanced, as explained below. 

Therefore, we use Ayps) = 0.039Aj/^( 15 ), A r ^ S ) = 0.039A T(1S) , and A r(3S ) = 
0.0064Ax(i5) in the standard model. On the other hand in the mixed hybrid 
study both ty'(2S) and T(35) were found to be approximately 50% hybrids. In 
Ref[8] it was shown, using the external field method, that the octet to singlet 
matrix element was enhanced by a factor of ir 2 compared to the standard model, 
as illustrated in Figure 3. For mixed hybrids we use an enhancement factor of 
3.0 




Figure 3: External field method for #'(25) and T(35) states 
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For differential cross sections we use the rapidity variable, y, 

1, / E + p z 



y (x) = — ^); with E = \J M 2 + p 



'z 



Pz = -rr(a:--), or (6) 



2 v x' 



x(y) = 0.5 



-^=(expy - exp (-y)) + J(^=(expy - exp (-y))) 2 + 4a 



(7) 



For the unpolarized proton collisions we use a polynomial fit to the parton 
distributions of Ref.[17]. Because of the wide range of vaues, in order to obtain 
a good polynomial fit to the parton distributions we limit the range of rapidity 

to -1. < y < I. 

For Q=3 GeV, with rn= Charmoiiium mass = 1.5 GeV, from Eq(7), x has a 
range about 0.028 to 0.032, and a/x 0.008 to 0.015. We have derived the following 
expressions for the gluon (g), u and d quark, and antiquark distribution functions 
using QTEQ6 for Q=3 GeV, fitting the range x=0.008 to .004, which is needed 
for v/i=200 GeV 



f g (x) ~ 1334.21 - 67056.5a; + 887962.0x 2 

f d (x) ~ 72.956 - 3281. Lr + 42247.6a; 2 

f u (x) ~ 82.33 - 3582.36a; + 45867.3a; 2 

f u (x) ~ 55.98 - 2722.04x + 35641. 2x 2 (8) 

f d (x) ~ 57.44 - 2757.05a; + 36030.5a; 2 . 

For Q=10 GeV, m=Bottomonium mass=5 GeV, from Eq(7), x has a range 
about 0.05 to 0.08, and a/x 0.03 to 0.05. We have derived the following ex- 
pressions for the gluon (g), u and d quark, and antiquark distribution functions 
using QTEQ6 for Q=10 GeV, fitting the range x=0.03 to .08, which is needed 
for v/i=38.8 GeV and 2.76 TeV. 

f g (x) ~ 275.14- 6167.6a; + 36871.3a; 2 
f d (x) ~ 26.96 - 527.14a; + 3119.13a; 2 
f u (x) ~ 32.92 - 604.38a; + 3530.1a; 2 

U(x) ~ 16.64 - 377.53a; + 2336.86a; 2 (9) 
f 3 (x) ~ 17.81 - 390.64a; + 2392.46a; 2 . 
The differential rapidity distribution for A = is given by 

= A^f g (x(y)^rn)f g (a/x(y),2m)^ , (10) 

while for A=l 

dapp ^ (x=1) = A*^\f a (x(y),2m)f g (a/x(y),2m) + 0MW 

dx 

+f u (x(y),2m)f a (a/x(y),2m)}— . (11) 

dy 
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3.1 Charmonium Production Via Unpolarized p-p Collisions at 
E =v ^= 200 Gev at BNL-RHIC 

First we consider unpolarized p-p collisions for y/s = 200GeV corresponding to 
BNL energy. We use scenerio 2 [15], with the nonperturbative matrix elements 
given above. Therefore, = §^f- g < Of^So) > =7.9 x l(T 4 nb for $=J/^ and 
2.13 x 10 5 nb for T(1S) heavy quark states; a = 4m 2 /s = 2.25 x 10~ 4 for Charmonium 
and 2.5 x 10 3 for Bottomium. 
For JTs = 200GeV 



x(y) = 0.5 
dx(y) M 



m 



(expy - exp (-y)) + 



dy 



400 



.200 

(expy + exp (-y)) 



(exp y — exp (— y))) 2 + 4a 



1. + 



m 
200 



||(expy + exp(-y)) 
^(exp y - exp (-y))) 2 + 4a 



(12) 



Using Eqs(10,ll,12), with the parton distribution functions given in Eq(8), 
we find da/dy for Q=3 GeV, A = and A = 1 the results for J/^> shown in Figure 
4. 



1.5 



1.0 



.5 





Figure 4: da/dy for Q=3 GeV, E=200 GeV unpolarized p-p collisions producing J/^fr with 
A = 0, A = 1 



Note that the shape of da/dy is consistent with the BNL-RHIC-PHENIX 
detector rapidity distribution[7]. 
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For ty'(2S) the results are shown in Figure 5. 



O.2. 




V|X ( 2S )(b)( A,=()) 
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V(2S)(a) (^=1) 

I 1 1 






— 1 
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Figure 5: da/dy for Q= 3 GeV, E=200 GeV unpolarized p-p collisions producing 1 !/' (2S) 
with A = 1,A = 



The results for dcr/dy shown in Figure 4. labeled ty'(2S)(a) are obtained by 
using for the standard nonperturbative matrix element =0.039 times the matrix 
elements for J/ty production; while the results labeled ty'(2S)(b) is obtained by 
using the matrix element derived using the result that the ty'(2S) is approximately 
50% a hybrid with the enhancement is at least a factor of n, as discussed above. 

3.2 Upsilon Production Via Unpolarized p-p Collisions at ~E=^/s= 
38.8 Gev at Fermilab 

Since the T(nS) states have not been resolved at BNL-RHIC at the present 
time, we consider T(nS) state production at 38.8, which has been measured at 
Fermilab[13, 14] 

For Q=10 GeV, using the parton distributions given in Eq(9) and Eqs(10,ll) 
for helicity A = 0, A = 1, with Ay =5.66 x lCT 4 nb and a = 6.64 x 1CT 2 , we obtain 
da/dy for T(nS) production. 
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The results for T(IS'), T(2S) are shown in Figure 6, and for T(3S) in Figure 7. 



Figure 6: du/dy for Q= 
T(2S) with A = 0, A = 



10 GeV, E=38.8 GeV unpolarized p-p collisions producing T(IS'), 



Figure 7: dcr/dy for Q= 10 GeV, E=38.8 GeV unpolarized p-p collisions producing T(3S) 
with A = 0, A = 1 
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In Figure 6 da/dy for T(1S) and T(2S) are obtained using the standard model 
for the matrix elements. In figure 7 the results for da/dy for the standard model 
are labelled T(3S')a, while the results for the hybrid model are labelled T(3S')b 

It should be noted that the ratios of da/dy for ty'(2S), shown in Figure 5 
and T(3S), shown in Figure 7 for the hybrid theory vs. the standard are our 
most significant results, as there are uncertainties in the absolute magnitudes 
and shapes of da/dy on the scenerios, as well as the magnitudes of the matrix 
elements. This is discussed in the following subsection. 

3.3 Dependence of da/dy on scenerios 

As stated above, our calculations make use of scenerio 2 of Ref[15]. It is impor- 
tant for us to point out that our main objective in the present work is to derive 
the relative magnitudes of the ty'(2S) and T(3S) cross sections in our mixed heavy 
hybrid theory, vs standard quark models for these states. To ilustrate this, we 
show Fig. 4 from Nayak and Smith' 2006 publication[15]. Scenerios 1, 2, 3 
are defined by 1 < Of( 3 P ) >= < Of (%) >= .0087GeV 3 , 2 < Of( 3 P ) >=0 and 
< OffSo) >=.039GeV 3 , and 3 < Of ( 3 P ) >= .01123Ge\/ 3 and < OfCS ) >= 0. See 
Ref[15] for the form of Eq(l) with the three scenerios. 




Figure 8: da/dy for Q= 3 GeV, E=200 GeV polarized p-p collisions with A = 1 producing 
J/ty. Scenerios 1, 2, 3 are shown with solid, dashed, and dotted curves, from Ref[15] 



As one can see from Figure 8, the shapes and magnitudes of da/dy depends on 
the scenerios, but this is not a problem in the present work, since we are mainly 
interested in ratios of cross sections; and our main problem with the shape (see 
Figs. 9 and 10) are for polarized ^ production as in Ref[15]. 
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4 Polarized p-p collisions at 200 Gev at BNL-RHIC 



For polarized p-p collisions the equations for °"'^ (A=U) and CTpp ^ (A=1) are the same 
as Eqs(10,ll) with the parton distribution functions fg and fq given in Eqs(8,9) 
replaced by Afg and Afq, the parton distribution functions for longitudinally 
polarized p-p collisions. 

A fit to the parton distribution functions for polarized p-p collisions for Q=3 
GeV obtained from CTEQ6[17] in the x range needed for ^=200 GeV is 

Af g (x) ~ 15.99 - 700.34a; + 13885.4x 2 - 97888.x 3 
Af d (x) ~ -5.378. + 205.60x - 4032.77x 2 + 28371.x 3 
Af u (x) ~ 8.44 - 292.19x + 5675.16x 2 - 39722.x 3 

Af a (x) ~ -1.447 + 64.67x - 1268.24x 2 + 8878.32x 3 (13) 
Af d (x) = Af u (x), 

and for Q=10 GeV, which we do not use in the present work, as the T(nS) are 
not resolved at BNL-RHIC, 



-6.074 + 241.57a; - 4762.04x 2 + 33604.4x 3 



A/ fl 10(x) ~ 28.98 - 1435.47x + 29533.5x 2 - 211440.x 3 
Af d l0(x 
A/ U 10(x 
A/ S 10(x 
A/jl0(x 



9.88 - 348.632x + 6729.49x 2 - 47058.x 



3 

-1.552 + 75.731x - 1531.97x 2 + 10896.6x 3 (14) 
A/ S 10(x) . 



The differential rapidity distribution for polarized p-p collisions are 

dAa pp ^ (x=0) = _^^ A/fl(x(y)j2m)A/ff(a/x(y)>2TO) ^ > (15) 



(IT 

Af 3 (a/x(y),2m) + Af u (x(y),2m)AU(a/x(y),2m)}— . (16) 

For polarized p-p collisions, Q=3 GeV, the results for dAcr/dy for J/ty pro- 
duction using the standard model are shown in Figure 9; while for ^>'(2S) the 
results are shown in Figure 10. As above, the curves labelled ^'(2S)a and are the 
standard model results, while that labelled ^J'(2S)b are the results for a mixed 
hybrid. The enhancement from active glue is once more quite evident. Since 
T(nS) states have not been resolved at BNL-RHIC, where polarized p-p collisions 
were measured, we do not calculate dAcr/dy for T(nS) states. 
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Figure 9: dAcx/dy for Q=3 GeV, E=200 GeV polarized p-p collisions producing J/\l/, with 
A = 0, A = 1 
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Figure 10: dAcr/dy for Q= 3 GeV, E=200 GeV polarized p-p collisions producing ty'(2S) 
with A = 0, A = 1 



Once again, we stress that it is the ratios of dAa/dy that is most significant, 
as there is uncertainty both the absolute magnitudes and shapes. Note that 
dAa/dy is expected to have a maximum at y=0. As in N-S for polarized J/ty 
production, dAa/dy < but there is a minimum at y=0 in absolute value for our 
results with scenerio 2. (See preceeding subsection on dependence on scenerios.) 
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5 Upsilon Production Via Unpolarized p-p Collisions 
at 2.76 Tev at LHC-CMS and 38.8 GeV at Fermilab 



The cross sections for T(nS) state production in p-p collisions have been mea- 
sured at 2.76 TeV at the LHC-CMS[12] and at 38.8 GeV at Fermilab [13]. In this 
subsection we calculate the cross sections for T(nS) production, with n= 1, 2, 
3. Then we use our theory that T(3S) is a hybrid to estimate the ratios of cross 
section. Since we are using scenerio 2 with < Of( 3 Po) >=0, the A = helicity 
dominates the cross section[15], and we drop the A = 1 terms. From Eq(5), for 
A = 0, the cross section is determined from 

o" PP ^$(A=o) = A* I — fg(x, 2m)fg(a/x, 2m) . 

J a X 

We use this to estimate the ratios of the T(2S) and T(3S) production cross 
sections to the T(IS') production cross section for 2.76 TeV and 38.8 GeV. 

5.1 Upsilon Production Via Unpolarized p-p Collisions at 2.76 
Tev 

For ^=2.76 TeV, for T(nS) production, 

a = 1.31 x 10" 5 
A r = 1.12 x 1(T 7 . (17) 

The cross sections for T(nS) state production in p-p collisions have been 
measured at 2.76 TeV at the LHC-CMS[12] and at 38.8 GeV at Fermilab [13]. 
In this subsection we calculate the cross sections for T(nS) production, with n= 
1, 2, 3. Then we use our theory that T(3S) is a hybrid to estimate the ratios 
of cross section. Since we are using scenerio 2 with < Of ( 3 Po) >=0, the A = 
helicity dominates the cross section[15], and we drop the A = 1 terms. From 
Eq(5), for A = 0, the cross section is determined from 



rl dx 

o- pp ^$(A=o) = Aq, — fg(x,2m)fg(a/x,2m). 

J a X 



rl dx 

,->-$(A=0) = Aq> I 
J < 

From this we find 



Opp^r(is) ^ 0.85 nb (18) 

What is significant for the present work are the ratios of the (IS), (2S), (3S) 
state production. In the preceeding sections we used estimates from the wave 
functions to find these ratios for the differential cross section. In the present 
section we make use of experimental results from the Fermilab experiment [13] 
for the ratio of cr(Y(2S'))/cr(Y(lS')), since the 2S and IS states are given by the 
standard model, and a(T(3S))/a(T(lS)) = N3 x a(T(2S))/a(T(lS)), as discussed in 
Section 2. Therefore, our estimate of the standard model is 
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a(T(2S))/a(T(lS)) ~ 0.27 standard 
a(T(3S))/a(T(lS)) ~ 0.04 standard, giving 
<j(T(2S)) + <r(T(3S)) 



<t(T(1S)) 



0.31 standard (19) 



On the other hand, in our mixed hybrid theory with the T(3S) about 50% 
hybrid [8], we would expect a factor of 7r 2 /4 in the matrix element, and therefore a 
factor of about 2.45 for the T(3S) cross section compared to the standard model. 

This results in our estimate 

*(?(2S)) + *(?(3S)) _ 

<t(T(1S)) " ' { } 

Compared to the LHC-CMS result [12] that this ratio is 0.78^ 6 14 ± .02. while in 
the standard model it would be about 0.31. 



5.2 Upsilon Production Via Unpolarized p-p Collisions at 38.8 
GeV 

Our study of the 38.8 Upsilon production is similar to the preseeding one for the 
LHC-CMS 2.76 TeV experiments. Our result for the a(T(3S))/a(T(lS)) expected 
at 38.8 GeV in the standard model, see Eq(19, compared to our mixed hybrid 
theory: 



gggg)) 

a(T(lS)) 
<r(T(3S)) 



0.04 standard 

0.147- 0.22 hybrid (21) 



a(T(lS)) 

compared to the experimental result[13] of about 0.12 to 0.16 



6 Conclusions 

We have applied the mixed hybrid theory for heavy quark states and predict 
that the cross sections for production of the charmonium ty'(2S) state in 200 
GeV p-p collisions and bottomonium T(3S) states in 38.8 GeV p-p collisions are 
much larger than the standard model. We have also estimated ratio of cross 
sections for 2.76 TeV and 38.8 GeV experiments, and our prediction for the 
1(3.5) production cross section is larger than the standard model, and closer to 
the experimental values. 

Because of the importance of gluonic production in processes in a Quark 
Gluon Plasma, this could lead to a test of the creation of QGP in RHIC. In order 
to treat the production of heavy mixed hybrid states via the QGP, however, we 
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must redo the QCD calculation using finite temperature field theory, as the 
properties of both the standard and hybrid components of the ty'(2S) and T(3S f ), 
as well as the properties of the J/ty and the other T(nS) states, are modified by 
the temperature during the QCD phase transition. This is a project for future 
research. 
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